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Abstract. In this paper we obtain a Hamilton-Jacobi theory for nonholo- 
nomic mechanical systems. The results are applied to a large class of nonholo- 
nomic mechanical systems, the so-called Caplygin systems. 



1. Introduction 

The standard formulation of the Hamilton-Jacobi problem for a hamiltonian 
system is look for a function S{t, q^) (called the principal function) such that 

f + H„^^, = 0. ,1.1, 

where H : T*Q — > R is the hamiltonian function. If one looks for solutions of the 
form S{t, q^) = VF(g^) — tE, where £' is a constant, then W must satisfy 

H{q\^)^E, (1.2) 

where W is called the characteristic function. 

Equations (jl.ip and p.2p are indistinctly referred as the Hamilton-Jacobi 
equation. 

The powerful of this method is that, in spite of the difhculties to solve a partial 
differential equation instead of an ordinary differential one, in many cases it works, 
being an extremely useful tool, usually more than Hamilton's equations. Indeed, 
in these cases the method provides an immediate way to integrate the equations of 
motion. The modern interpretation relating the Hamilton-Jacobi procedure with 
the theory of lagrangian submanifolds is an important source of new results and 
insights [2 [2]. Let us remark that, recently, Carihena et al [5] have developed a 
new approach to the geometric Hamilton-Jacobi theory. 

On the other hand, in the last fifteen years there has been a renewed interest in 
nonholonomic mechanics, that is, those mechanical systems given by a lagrangian 
L = L{q^,q^) subject to constraints ^^{q^,q^) = involving the velocities (see [5] 
and references therein). A relevant difference with the unconstrained mechanical 
systems is that a nonholonomic system is not hamiltonian in the sense that the 
phase space is just the constraint submanifold and not the cotangent bundle of 
the configuration manifold; moreover, its dynamics is given by an almost Poisson 
bracket, that is, a bracket not satisfying the Jacobi identity [S]. In [TT] the au- 
thors proved that the nonholonomic dynamics can be obtained by projecting the 
unconstrained dynamics; this will be the point of view adopted in the present paper. 
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A natural question, related with a possible notion of integrability is in what 
extent one could construct a Hamilton- Jacobi theory for nonholonomic mechanics. 
Past attempts to obtain a Hamilton-Jacoby theory for nonholonomic systems were 
non-effective or very restrictive (see [71 [T9j [201 Ell l22^ and also ^TSj), because, in 
many of them, they try to adapt the typical proof of the Hamilton- Jacobi equations 
for systems without constraints (using Hamilton's principle). Usually the results 
are valid when the solutions of the nonholonomic problem are also the solutions of 
the corresponding constrained variational problem (see |101 117| I18j for a complete 
discussion) . 

In our paper, we present an alternative approach based on the geometrical prop- 
erties of nonholonomic systems (see also [16] for second-order differential equa- 
tions). The method is applied to a particular class of nonholonomic systems, called 
Caplygin systems: in such a system the configuration manifold is a fibration over 
another manifold, and the constraints are given by the horizontal subspaces of a 
connection on the fibration. In this case, the original nonholonomic system is equiv- 
alent to another one whose configuration manifold is the base of the fibration and, 
in addition, it is subject to an external force [H]. In any case, the equations we 
obtained are different that in previous works and may give new insight in this topic. 
In particular, this theory could give insights in the study of integrability for non- 
holonomic systems [3] and even in the construction of new geometrical integrators 
for nonholonomic systems (see [UlfT^). 

2. Preliminaries 

2.1. Lagrangian and Hamiltonian mechanics. Let L = L{q'^,q'^) be a la- 
grangian function, where {q^) are coordinates in a configuration n-manifold Q. 
Hamilton's principle produces the Euler-Lagrange equations 

d f dL\ dL , , , 

A geometric version of Eq. (|2.ip (see [T^) can be obtained as follows. Consider 
the (l,l)-tensor field S and the Liouville vector field A locally defined on the tangent 
bundle TQ of Q by 

d , A ^ -A 9 

oq^ oq^ 
Since the lagrangian L is a function defined on TQ one can construct the Poincare- 
Cartan 1- and 2-forms 

aL = S*{dL), ujl = -daL, 

where S* denotes the adjoint operator of S. The energy is given by El — A(L) —L. 
We say that L is regular if the 2-form uj^ is symplectic. In this case, the equation 

ix UJL = cLEl (2.2) 

has a unique solution, X — called the Euler-Lagrange vector field; is a second 
order differential equation (SODE) which means that its integral curves are tangent 
lifts of their projections on Q (these projections are called the solutions of ^l). A 
direct computation shows that the solutions of are just the ones of Eqs. (|2.ip . 

Finally, let us recall that the Legendre transformation FL : TQ — > T*Q is a 
fibred mapping (that is, ttq o FL = tq, where tq : TQ — > Q and ttq : T*Q — > 
Q denote the canonical projections of the tangent and cotangent bundles of Q, 
respectively). The regularity of L is equivalent to FL being a local diffeomorphism. 
Along this paper, we will assume that FL is in fact a global diffeomorphism (in other 
words, L is hyperregular) which is the case when L is a lagrangian of mechanical 
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type, say L — T ~V, where T is the kinetic energy defined by a Riemannian metric 
on Q and V : Q — > R is a potential energy. 

The hamiltonian counterpart is developed in the cotangent bundle T*Q of Q. 
Denote by luq = dq^ A dpA the canonical symplectic form, where are the 

canonical coordinates on T*Q. The Hamiltonian energy is just H = El o FL~^ 
and the Hamiltonian vector field is the solution of the symplectic equation 

As we know, the integral curves [q^ [t) , p A{t)) of Xh satisfy the Hamilton equations 

dH 



OPA 

dH 



(2.3) 



Finally, since FL*loq — ujl we deduce that and Xh are _FL-related and, 
consequently, FL transforms the Euler-Lagrange equations (|2.ip into the Hamilton 
equations 



2.2. Nonholonomic mechanical systems. A nonholonomic mechanical system 
is given by a lagrangian function L = L{q^,q"^) subject to a family of constraint 
functions 

<i>'{q^,q'^) =0, 1 < i < TO < n = dimQ. 

In the sequel, we will assume that the constraints are linear in the velocities, 
i.e., $^(q^,q^) = <P\iq)q^. 

Invoking the D'Alembert principle, we derive the nonholonomic equations of 
motion 

(p^^qA^A) = 0, l<i<TO J 

where Xi ~ Xi{q'^, g'^), 1 < i < m, are Lagrange multipliers to be determined. 

In a geometrical setting, L is a function on TQ and the constraints are given by 
a vector subbundle M of TQ locally defined by = 0. 

Equations (|2.4p can be intrinsically (see jT5]) rewritten as follows 



ixLOL-dEL e S*{{TMf) 
X e TM. 



(2.5) 



For the formulation of a Hamilton- Jacobi theory we are interested in the "Hamil- 
tonian version" of the nonholonomic equations. Assuming that the Lagrangian L 
is hyperregular, then the constraint functions on T*Q become = o FL~^, i.e. 

r) H 

*Xg-',PA) = <i>k(9)^(9^,PA), 

OpA 

where the Hamiltonian H : T*Q is defined hy H — El o FL^^. 

The equations of motion for the nonholonomic system on T*Q can now be written 
as follows 

.A ^ dH_ >| 

- - > (2-6) 



PA = - \.<i>\{q) 

together with the constraint equations ^''(q,^) = 0. 
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Let M denote the image of the constraint submanifold M under the Legendre 
transformation, and let F be the distribution on T*Q along M, whose annihilator 
is given by 

F° ^ FL,{S*{{TMf)). 

Observe that F'^ is locally generated by the m independent 1-forms 

^« = ^\{q)dq^ , 1 <i <m. 

The nonholonomic Hamilton equations for the nonholonomic system can be then 
rewritten in intrinsic form as 



\M 

e TM 



(2.7) 



Assume the compatibility condition: F^ n TM = {0}, where " _L " denotes 
the symplectic orthogonal with respect to loq. Observe that, locally, this condition 
means that the matrix 

(e^) = (<i>^(g)J{^^<i>i3(g)) (2.8) 

is regular, where {'K'^^) = {d^H/dpAdps)- The compatibility condition is not 
too restrictive, since it is trivially verified by the usual systems of mechanical type 
(i.e. with a Lagrangian of the form kinetic minus potential energy). The compat- 
ibility condition guarantees in particular the existence of a unique solution of the 
constrained equations of motion (j2.7p which, henceforth, will be denoted by Xnu- 

Moreover, if we denote by Xh the Hamiltonian vector field of if, i.e. ix^^Q = 
dH then, using the constraint functions, we may explicitly determine the Lagrange 
multipliers as 

A, = eyXH(*^) , (2.9) 
where (Cy ) is the inverse matrix of (C*^). 

2.3. Caplygin systems. A Caplygin system is a nonholonomic mechanical system 
such that: 

(i) the configuration manifold Q is a fibred manifold, say p : Q — > N, over a 
manifold N; 

(ii) the constraints are provided by the horizontal distribution of an Ehres- 
mann connection T in p; 

(iii) the lagrangian L : TQ — > R is F-invariant. 

Remark 2.1. A particular case is when p : Q — > N = Q/G is a principal G-bundle 
and F a principal connection. o 

Let us recall that the connection F induces a Whitney decomposition TQ = 'K(BVp 
where IK is the horizontal distribution, and Vp = ker Tp is the vertical distribution. 
Take fibred coordinates (g'^) = {q'^,q^) such that p{q'^,q^) — (g"); therefore we can 
obtain an adapted local basis of vector fields 

^= {^a = 1^ - K^), Vp={V, = 4-)- 
^ dq" "dq"' ^ ^ dq"-' 

d d 

Here 5{a = (— — )^ = h(— — ), where denotes the horizontal lift of a tangent 

vector y on iV to Q, and h : TQ — > % is the horizontal projector; F^ — T\{q^) 
are the Christoffel components of the connection F. 

The dual local basis of 1-forms is 

{ila^dq\'n,=dq'+Tldq''} 
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The curvature of F is the (l,2)-tensor field R= ^[h, h] where [ , ] is the Nijenhuis 
tensor of h, that is 

R{X, Y) = [hX, hY] - h[hX, Y] - h[X, hY] + h^[X, Y] 

Therefore we have 



where 



d d d 



^ab Q^b Qqa + a q^j b g^j 



The constraints are locally given by (f>* = (7* + T^g" = 0. In other words, the 
solutions are horizontal curves with respect to F. 

Since the lagrangian L is F-invariant, that is, L{(Y^)q-^) = L{{Y^)q^) for all 
Y £ TyN, y = p{qi) = p{q2), we can define a function L* : TN — > M as follows: 
L*0^v) — L{0^^)q)^ where y = p{q). Therefore we have 

L*{q-,r)^Liq%q\r,'~Kr) 



(2.10) 



Equations (|2.5p read now as 

ixuJL'dEL e 5*((T5{)") 
X (E TJ{ 

Define a 1-form a* on TN by putting 

where U e r„(riV), u e TyN, U G T^iTQ) such that C7 projects onto 

R{{u)^,{TrN{U)^))eTqQ, 
p{q) = y, X G Tgi^) = Q- In local coordinates we obtain 

^'i'Kb) dqr 



Consider the following equation 

iyujL' - dEL' = a*. (2.11) 

A long but straightforward proof shows that L* is a regular lagrangian on TN , 
therefore (|2.1ip has a unique solution Y* . Notice that the pair (L*,a*) can be 
considered as an unconstrained system subject to an external force a* . The corre- 
sponding equations of motion are 

fKb- (2.12) 

Both systems, the nonholonomic one on Q given by L and the constraints given 
by F, and that given by L* and a*, are equivalent. The equivalence is explained in 
the following. 

F induces a connection F in the fibred manifold Tp : TQ — > TN along !K by 
defining its horizontal distribution as follows: 



d fdL*\ 


dL* 


dL 


dt V J 


dq°- 


dqi 



f—)-^ ^ — p — - a — ^ - r-'— i — 

^Qqa^ Qqa a Q^i y Q^a a g^, ) 
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Theorem 2.2. The nonholonomic dynamics Xnh is a vector field on 5f which is 
Tp-projectable onto Y*. Furthermore, Xnh is the horizontal lift ofY* with respect 
to the induced connection T. 

Example 2.3. (Mobile robot with fixed orientation) 

The body of the robot maintains a fixed orientation with respect to the environ- 
ment. The robot has three wheels with radius R, which turn simultaneously about 
independent axes, and perform a rolling without sliding over a horizontal floor. 

Let {x, y) denotes the position of the centre of mass, 6 the steering angle of the 
wheel, "0 the rotation angle of the wheels in their rolling motion over the floor. So, 
the configuration manifold is Q = S*^ x 5^ x K^. The lagrangian L is 

1 1 1 • 3 • 

L = -mx^ + -mf + -J9^ + -J^iP^ 

where m is the mass, J is the moment of inertia and J^^ is the axial moment of 
inertia of the robot. 

The constraints are induced by the conditions that the wheels roll without sliding, 
in the direction in which they point, and that the instantaneous contact point of 
the wheels with the floor have no velocity component orthogonal to that direction: 

i: sin — y cos = 0, 
X cos 9 + y sin 9 — Rip — 0. 

The abelian group G = acts on Q by translations, say 

((a, b), (9, ip, X, y)) <^ (9,'i\),a^ x,b^ y) 

Therefore we have a principal G-bundle p : Q — > N ~ Q/G with a principal 
connection given by the connection 1-form 

P — (dx — Rcos9d'i(j)ei + [dy — Rfi\n9dxp)e2 

where {ei, 62} denotes the standard basis of R^. The constraints are given by the 
horizontal subspaces of 13. If we apply the above reduction procedure we deduce 
a* = 0. 

3. Geometric Hamilton- Jacobi theory 

The following result is a geometric version of the standard formulation of the 
Hamilton- Jacobi problem [I]. 

Theorem 3.1. Let j be a closed 1-form on Q. Then the following conditions are 
equivalent: 

(i) for every curve a : R — > Q such that 

&{t)^TnQ{XH{j{<Jm) 

for all t, then ^ o a is an integral curve of Xh- 

(ii) d{Ho-f) = 0. 

If 7 = dW we recover the standard formulation since d{HodW) = is equivalent 
to the condition H o dW = cte, that is 

where E is a. constant. 

A interesting new point of view of the geometric Hamilton- Jacobi theory has 
been recently developed by J.F. Carihcna et al. 
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Let 7 be a closed 1-form as in Theorem 13.11 Since FL is a diffeomorphism, we 
can define a vector field X on Q by 

X = FL^i 07 

Therefore, we have 

= d{H 07) = d{EL o FL-^ 07) = d{EL o X) 
because H = E^o FL^^. 

Hence, Theorem 13 . 1 1 can be reformulated as follows. 

Theorem 3.2. [B] Let X he a vector field on Q such that FLoX is a closed 1-form. 
Then the following conditions are equivalent: 

(i) for every curve a : M. — > Q such that 

<y{t)^TrQ{UX{am) 

for all t, then X o a is an integral curve of 

(ii) d{ELoX) = 0. 

Definition 3.3. A vector field X satisfying the conditions of Theorem \3.2\ will be 
called a solution for the Hamilton- Jacobi problem given by L. 

3.1. An interlude: meclianical systems witii external forces. We shall need 
the following formulation of the Hamilton- Jacobi theory for mechanical systems 
with external forces. 

A mechanical system with an external force is given by (see [S]): 

(i) A lagrangian function L : TQ — > M, where Q is the configuration mani- 
fold; 

(ii) a semibasic 1-form a on TQ. 

Since a is semibasic (that means that a vanishes when it is applied to vertical 
tangent vectors) we have 

a = ctAiq, q)dq^ 

The Euler-Lagrange equations are then 
d dL dL 
dt^d^' ^ dq" 
which correspond to the symplectic equation 

ix iOL = dEL+a (3.2) 

Indeed, when L is regular, Eq. p.2p has a unique solution ^^^q, which is a second 
order differential equation whose solutions are just the ones of p.ip . 

Notice that ^l,q = ^.l + Z, where iz^L = ot- 

Observe that we can construct the hamiltonian counterpart using the Legendre 
transformation, so that we have a hamiltonian H = El ° FL~^ subject to the 
external force /3 = {FL^^)*a which is again semibasic (i.e. (3 — (3Adq^). The 
equation 

^^Q^dH + 13 

has a unique solution Xn.p whose integral curves satisfy the Hamilton equations 
with external force 

■A 9H •) 

Q = — I 

(3.3) 



-iinA)-7^ = -^A,l<A<n, (3.1) 



PA 
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Theorem 3.4. Let ^ be a closed 1-form on Q. Then the following conditions are 
equivalent: 

(i) for every curve a : M — > Q such that 

&{t)=TnQ{XHAlim})) (3.4) 

for all t, then "f o a is an integral curve of Xh,!)- 

(ii) d{H 07)== -7*/3. 

Proof. Since 7 = 7_a dq"^ is closed then 

d"fA _ djB 
dq^ dq^ 

It is easy to show that Equation p.4p is rewritten, in local coordinates, as 



&^{t) = j—{<J^{t),^B{am (3.5) 

OPA 



We also have that condition 

07) = -7*/3 

is written in local coordinates as 

dH dH djB 



dq^ dps dq"^ 
(=4>) Assume that (i) holds. Therefore 



- -f^A (3.6) 



J n rr 

-(7A(a(t))) = -_(7(a(0)) - PAhHt))) (3.7) 



Moreover 



dH dH djB dH dH d^A , . . , ,^ 

— ' (since 7 IS closed) 



dq^^ dpB dq^ dq^ dpB dq^ 

= -/3a (from (1X71) 1 
(<^=) Assume that (ii) holds, that is. 



dH dH djB 
dq^ dpB dq^ 
Now using p.Sp and since 7 is closed, then 



A 



Therefore 



dH ..Bf.^dlA 



J^ilAHt))) = -^{^{am-PAilHt))) 
which proves that 7 o cr is an integral curve of Xu.p- D 

Therefore we have the lagrangian version. 

Theorem 3.5. Let X be a vector field on Q such that FL o X is a closed 1-form. 
Then the following conditions are equivalent: 

(i) for every curve a : M — > Q such that 

for all t, then X o a is an integral curve of ^L,a- 
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(ii) d{ELoX) = -X*a. 

Definition 3.6. A vector field X satisfying the conditions of Theorem \3.5\ will be 
called a solution for the Hamilton- J acohi problem given by L and a. 

4. Hamilton- Jacobi theory for nonholonomic mechanical systems 

Let L : TQ — > M be a lagrangian function subject to nonholonomic constraints 
given by a vector subbundle M of TQ, locally defined by the linear constraints 

— 1 < i < m. Denote by D the distribution on Q whose annihilator 

is L>° = span{^* = <l>\{q)dq'^}. Notice that S*{TM°) is the puUback to TQ of the 
annihilator _D" of D. 

We assume the admissibility and compatibility conditions, and consider the 
hamiltonian counterpart given by a Hamiltonian function H : T*Q — > R and 
a constraint submanifold M = FL{M) as in the precedent sections. Xnh and Xnh 
will denote the corresponding nonholonomic dynamics. Given D^, the annihilator 
of D, we can form the algebraic ideal 'J{D'~') in the algebra A*{Q). Therefore, if a 
fc-form 1^ e 3{D°) then 

ly = Pi A ^\ where Pi £ A''"^((5), I < i < m. 

Theorem 4.1. Let j be a 1-form on Q such that j{Q) C M and dj £ Then 
the following conditions are equivalent: 

(i) for every curve a : R — > Q such that 

a{t)^T7TQ{XHh{<T{tm} (4.1) 

for all t, then "f o a is an integral curve of Xnh- 

(ii) d(Hoj) e D". 

Proof. The condition dj £ 3(1?°) means that 

where 7 = '^Adq^ and Pi — PiAdq^. It is easy to show that Equation (|4.ip is 
rewritten, in local coordinates, as 

pi T7 

crHt)^—{a^{t),^B{<j{t))) (4.2) 

OPA 

We also have that condition 

d{H o 7) e D'^ 

is written in local coordinates as 

■ dH dH djB 



dq"" = K^\{q)dq'' (4.3) 



dq"^ dpB dq 

for some Lagrange multipliers A^'s. 

(=>) Assume that (i) holds. Therefore 

d f)M - 

-(7A(a(t))) = -_(7(a(t))) - X.^\{a{t)) , (4.4) 
where the Aj's are determined using the constraint equations 

^\a{t),^{am = ia{t),j{a{t)WBHt)) = 0. 

OPB 
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Using the constraint equations we deduce that 

dH dH d-fB dH dH djA OH _ dH 



dq"^ ^ dpB dq-^ dq^ ^ dps dq^ ^ dps^^^'^^ dp"^^^^^ 

dH 



dpi 

Therefore, we conclude that d{H o 7) G 
(<^=) Assume that (ii) holds, that is, 

'dH dH d-fB 
dq^ dpB dq^ 

Now using (|4.ip and since dj e 'J{D'^), then 



P^B (from 



dq^ = ~XiH' 



A 



l{q) = [q,— ~~Kil\), (4.6) 



Therefore 

|(7a(^(0)) = -^(7(^(0)) + (a. - a^mBHt))) ^\{a{t)) (4.5) 

Using that Im(7) e M, we deduce that \i = ■§^PiB — K along 7. □ 

Remark 4.2. Suppose that j — dS where 5 is a function S : Q — > M. In this case, 
the condition d'y G 3{D^) is trivially satisfied. Moreover, we note that in previous 
approximations to Hamilton- Jacobi theory [71 [TH [TH [171 [18] the considered sections 
are of the form 

dS 

and the coefficients Xi are determined through the nonholonomic constraint equa- 
tions 

dH 

Ma(9)^('7,7a(«)) = . 

In general, this type of 1-forms does not satisfy the condition that we initially 
impose, d'j G 3{D'^). Observe that in the particular case of holonomic constraints 
both approaches coincide. o 

Now, we write a coordinate expression for the Hamilton- Jacobi equation that 
we have proposed. In order to do it, consider a set of independent vector fields 
{Za = Z^-^}, 1 <a < n-m, on Q such that p''{Za) = 0, i.e, Dg = span {{Za)\q}. 
Thus a 1-form 7 on Q, solution of the nonholonomic Hamilton- Jacobi equation, 
must verify the condition d'f £ 3{D°) and, additionally, 

^.^('?)(^('?'7('/)) + ^(9,7(<z))f^(9)) = 0, 

P^'A{q)g^{q,l{q)) = 0, 

for the condition d{H o 7) G 1?° and for the condition jiQ) C M, correspondingly. 

Theorem 4.3. Let X he vector field on Q such that X{Q) C M and d{FL o X) G 
^{D^). Then the following conditions are equivalent: 
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(i) for every curve a : R — > Q such that 

a{t)^TTQ{X,,h{X{a(t)))) (4-7) 

for all t, then X o a is an integral curve of Xnh- 

(ii) d{ELo X) e 

Definition 4.4. A vector field X satisfying the conditions of Theorem \4.1\ will be 
called a solution for the Hamilton- Jacobi problem given by L and M . 

4.1. An application to Caplygin systems. Consider now the case of a Caplygin 
system (see Section 12. 3p . That is, we have a fibration p : Q — > iV, and an 
Ehresmann connection F in p, whose horizontal distribution imposes the constraints 
to a lagrangian L : TQ — > R. 

Let L* : TN — > M be the reduced lagrangian and a* the corresponding external 
force. We denote by Xnh the nonholonomic vector field on TQ and by X* the 
solution of the reduced lagrangian system with external force a* . 

Theorem 4.5. (i) Assume that a vector field X on Q is a solution for the 

Hamilton- Jacobi problem given by L and T. If X is p-projectable to a 
vector field Y on N and d{FL* o Y) — then Y is a solution of the 
Hamilton- Jacobi problem given by L* and a* . 
(ii) Conversely, let Y be a vector field which is a solution of the Hamilton- 
Jacobi problem given by L* and a*. Then, if d{FL o Y^) G 3(!K°), the 
horizontal lift Y^ is a solution for the Hamilton- Jacobi problem given by 
L and p. 

Proof (=^) 

Assume that a vector field X on Q is a solution for the Hamilton- Jacobi problem 
given by L and F, and that X is p-projectable onto a vector field Y on N. We have 
to prove that Y is then a solution of the Hamilton problem given L* and a* . Let 
p a curve in N such that 

m=Tm{Y*iYipitm (4.8) 

for all t. Take an horizontal lift a oi p to Q with respect to the connection F. A 
direct computation shows that 

a{t)=TTQiXnh{X{am) (4.9) 

since Xnh is the horizontal lift of Y* with respect to the prolongated connection F. 
Therefore we have that A o (j is an integral curve of Xnh and, consequently, Y o p 
is an integral curve oiY*. 

(^) 

Assume that Y is vector field on N which is a solution of the Hamilton- Jacobi 
problem given by L* and a*. Take its horizontal lift X = Y^ to Q. with respect 
to F. If cr is a curve in Q satisfying 

a{t)^TTQiXnh{X{am) (4.10) 

then the projection p = p o a satisfies ()4.8p . So, Y o p \s an integral curve of Y* 
and, hence A o cr is an integral curve of Xnh ■ 

□ 

Example 4.6. (The mobile robot with fixed orientation revisited) The reduced 
lagrangian in this case is 

L {0,'ip)^-J9 H ?A , 
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and 



and a* = 
Therefore, 

arc solutions of the Hamilton- Jacobi problem given by {L*,a*). Calculating the 
horizontal lifts off both vector fields wc have that: 

Y-^ = ^ and y« = A + Rcose^ + RsinO^ . 
of) dtp ox oy 

Now 

71 = FLoY^ = Jde 

72 = FL o Y^ = 3 Jo, dtp + mR cos 6dx + mR sin 6 dy 

^72 = -mRde A (sin Odx- cos 6 dy) e 

Therefore, Y^ and are solutions of the Hamilton-Jacobi problem of the non- 
holonomic problem given by {L, !K). Observe that in both cases d{H o 7^) = 0, for 
i = 1,2. In such a case, 

t I — > (xo,yo,t + eo,tpo) 

t I — > {tRcosOi} + XQ,tRsm6a + yi},6i:nt + tpa) 

are the solutions of the nonholonomic system {L,'K) obtained from Y^ and Y^ , 
respectively. Both solutions are also solutions of the lagrangian system determined 
by L without constraints; indeed, they are solutions of the free system satisfying 
additionally the nonholonomic constraints. 

But taking now the vector field 

it is obviously a solution of the Hamilton-Jacobi equations for the lagrangian L* 
and its horizontal lift 

is a solution of the Hamilton-Jacobi equations for the nonholonomic system {L,H): 

73 = FLoY^ = Jd6 + 3J^dtp + mRcos0dx + mRsm0dy 

and ^73 € In such a case, the solution of the nonholonomic problem that 

we obtain is 

t I — > {Rsin{t — 9o) + xq + Rsm9o,~Rcoa{t — 9o) + Vo + Rcos9o,t + 9o,t + tpo) 
which is a solution of the nonholonomic problem but not of the free system. 
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